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The interelectronic-interaction effect on the transition probabilities in high-Z He-like ions is in- 
vestigated within a systematic quantum electrodynamic approach. The calculation formulas for the 
interelectronic-interaction corrections of first order in 1/Z are derived using the two-time Green 
function method. These formulas are employed for numerical evaluations of the magnetic transition 
probabilities in heliumlike ions. The results of the calculations are compared with experimental 
values and previous calculations. 

PACS numbers: 32.70.Cs, 31.30.Jv 

I. INTRODUCTION 

During the last few years, transition probabilities in heliumlike ions were calculated by a number of authors P, 0, 
1^^. In these calculations, the interelectronic-interaction effects on the transition probabilities were accounted for 
by employing the relativistic many-body perturbation theory (RMBPT) Q, |^ Q and the multiconfiguration Dirac- 
Fock (MCDF) method pj. Since these methods are based on using the Coulomb-Breit hamiltonian, they have to 
deal with a separate treatment of the positive- and negative-energy state contributions. As was first indicated in 
0, the contribution from the negative-continuum contribution is very sensitive to the choice of the one-electron 
model potential, which is used as the starting point of any RMBPT or MCDF calculation. In particular, using a 
standard Dirac-Fock approximation, in Q it has been demonstrated that to achieve the agreement between theory 
and experiment for the magnetic dipole transition 2^ Si — > l^S'o it is necessary to take into account both correlation 
and negative-continuum effects. This statement is closely related to a problem of significant numerical cancellations 
that may occur in low-Z systems, if an improper one-electron approximation is used. For a rigorous QED approach 
for low-Z systems and corresponding calculations we refer to 

The main goal of the present paper is to perform a complete QED calculation of the interelectronic-interaction 
correction of first order in 1/Z to the magnetic transition probabilities in high-Z He- like ions. To derive the calculation 
formulas for these corrections from the first principles of QED we use the two-time Green function method developed 
in 0,13,13 and described in details in 0. In Sec. ^ we formulate the basic equations of this method for the case of 
nondegenerate states and apply it for the derivation of the desired formulas. The numerical results for the transitions 
2^Si 1^5*0, 2^P2 l^S'o, and 3'^ 6*1 2^Si are presented in Sec. IIIII Both Feynman and Coulomb gauges are used 
for the photon propagator to demonstrate the gauge independence of the final results. The results of the calculations 
are compared with previous theoretical results and with experiment. 

The relativistic units {h = c = 1) and the Heaviside charge unit ( a — e^/(47r), e < ) are used in the paper. 

II. BASIC FORMULAS 

We consider the transition of a high-Z two-electron ion from an initial state a to a final state b with the emission 
of a photon with momentum k / and polarization e f . The transition probability is given as 

dW = 2TT\T^f^b;af6{Eb + k"} - Ea)dkf , (1) 

where T^f^b-a is the transition amplitude which is connected with the 5- matrix element by 

(k/, e/; 6|5|a) = 27riT.^j.^b;a'5(£6 + - £a) , (2) 

and fcj. = |k/|. 



We assume that in zeroth (one-electron) approximation the initial and final states of the ion are described by 
one-determinant wave functions 



ia(xi,X2) = --^^(-l)^'(/'Pai(xi)V'Pa2(x2), 



Ub 



(X1,X2) = -^^(-l)'^VP6i(xi)^Pb2(x2). 



(3) 
(4) 



To describe the process under consideration we introduce the Green function gy^^b;a{E' , E) by 
g^,,bAE', E)5{E' + k'-E) = - J dpldpldp'^ dp'^ 

x5{E~pI-pI)5{E' ~p'^-p'^) 
X j (ixidx2dx']^dx2u|,(x']^, X2) 

X G,^ ( (pf , x; ) , ip'^ , x^ ) ; fcO ; (pO , xi ) , , X2 Mf^^a (xi , X2 ) , (5) 



where 



G^,((p'°,x;), (p^°,x^); fcO; (p;,Xi), {pl^^)) 



2-K 1 
T (27r)5 



X exp iip[^xf + ip'^^x'^' - ip'ixl - ip'ix'i -f ik^y'^) 
X Ay (y) (0|TV'(x'i (y)V^(x2)V^(xi ) |0) 



(6) 



is the Fourier transform of the four-time Green function describing the process, ■(/'(x) is the electron-positron field 
operator in the Heisenberg representation, and 



AKx) 



exp [ikf ■ x) 



(7) 



/2fc0(27r)3 

is the wave function of the emitted photon. The transition amplitude S^j^b;a = (k/, e^; 6|5|a) is calculated by 0,0,01 



S. 



7/, 



Z. 



-1/2 



6{Eb + kj ~ Ea) dE' j dEg^f^b.a{E',E) 



r 1 



dEgbbiE) 



-1/2 



1 

2Td 



dEgaa{E) 



'1/2 



(8) 



Here gaa{E) is defined by 
gaa{E)5{E' - E) = 



dp\dpldp'^dp'^ 



2tt 1 
T2! 

X^(i?-p?-p^)5(i?'-pf-p^0) 

X / dxidx2dx;dx^ x4(x'i,x^)G((p'°,x;),(pi°,x^);(p;,Xi),(pO,X2))7?72V(xi,X2), (9) 



where 



G(bf,x;),(p^°,x'2);(p?,xi),(p°,X2)) 

^ / da:'2° exp {ipfxf + ip'°x'° - ipl 



(27 

x(0|TVK)V(4)V'(a:2)V'(^i)|0) 



Xi IP2X2) 



(10) 



is the Fourier transform of the four-time Green function describing the ion; gbb{E) is defined by a similar equation. 
The contours Fa and Vb surround the poles corresponding to the initial and final levels and keep outside all other 
singularities of the Green functions. It is assumed that they are oriented anticlockwise. The Green functions G and 
G^fj are constructed by perturbation theory after the transition to the interaction representation and using Wick's 
theorem. The Feynman rules for G and G^^ are given, e.g., in 10]. 

Below we consider the transition probability in high-Z He- like ion to zeroth and first order in 1/Z. 



A. Zeroth order approximation 

To zeroth order in XjZ the transition amphtude is described by the diagrams shown in Fig. 1^ Formula (jS]) gives 

Cm = ^^^^ + ^) / / ^) ' (11) 

where the superscript indicates the order in XjZ. According to definition jSJ and the Feynman rules for G-y^ |lC| . we 
have 

/OO 



X 



« V- I»2)(n2| I , , I i l'^3)(»3| I m On 



"1 



"3 



where a'^ = -^O-y^ — (!,«)■ One obtains 



- 2 /" c'o -1 1 A 

5^°) (£;',S) = f— W (fe^lea^^flai)^^ ^ 



1 



where -E^"^ = + £02 ^-nd i?^"'' = e^^ + ■ Substituting this expression into equation and integrating over E 
and E' we find 

C6;a = -27rz(5(Sfc + fc° - E,)[(6i|ea^^f |ai),5,,fc, + (^zjea^^f |a2)5aifcj (15) 
or, according to definition ||2Jl, 



,(0) 

The corresponding transition probability is 



'^^I'.b-.a = -[{bi\ea^,A^j*\ai)5a,_b2 + (^2|ea^A^*|a2)(SaibJ . (16) 



B. Interelectronic-interaction corrections of first order in 1/Z 

The interelectronic-interaction corrections to the transition amplitude of first order in 1 / Z are defined by diagrams 
shown in Fig. 2a, b. Formula ^ yields in the order under consideration 

= 5{E, + k°f-Ea) U dE' f dEg^l\jE',E) 

'""'i '^^Ca^'^^) {hi dEg^^iE))] , (18) 



where gi]}{E) and {E) are defined by the first order interelectronic-interaction diagram (Fig. 3). Let us consider 
first the contribution of the diagrams shown in Fig. 2a. According to the definition @ and the Feynman rules for 
G-if we have 

g':^;lSE\E)5{E' + k^~E) 



/oc 
dpldpldp'^dp'^5{E -pI- pI)5{E' pf p'^ 
j_. -oo 

■ \ 3 f 1 



1 



'gO-£„(l-iO) 

x(nP6.|/(.)|a,a.)^„ - + zOp- - U + ^0 p", - e\ + ^0 

x<5(rf -oj- (7°)5(<z" ~k^~ )<5bO + c - p^O) 

+ -;n y^(Pb2\eanA't*\n)— \ ^ 

x(P&in|/(u;)|aia2)— 111 



p§ - £a2 + *0 - £pb, + iO - + iO 

x5{p"^ -Lo- q'')S{q° - fcO - + - )} , (19) 

where I{uj) = e^Q;^a''-Dpi/(ti;) and 

^ , ^ f d]i exp (ik • (x — y)) 

D,.i., X - y) . y ^ ^Li.. + ,o ^^^^ 

is the photon propagator in the Feynman gauge. One finds 

1 1 
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' pl-Ea^ +iOE-p^-ea, +iO 
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^ -{nPb2\I{p'2'-p'2)\aia2) 



E - pf - en{l - tO) 



p „ .-oo p'l" - £Pbi + io - p'l" - £pf„ + io 

^ ^ (P&2|ea^Af |n) 



^ (P6i7i|/(pf-p?)|aia2). (21) 



^-p'0-£„(l-zO) 

The expression (|21|l is conveniently divided into irreducible and reducible parts. The reducible part is the one with 
£pb2 + £n = in first term and with Ep^j^ + £n = E^'^ in second term. The irredicible part is the reminder. Using 
the identities 

~W) («?-£., +zO + S -«?-£„„ +zo) ' ^^^^ 



p? -£ai +iOS-pO -£a2 +iO £; - £;(°npO - £<,^ + iO S-pO-Ea^+iO/ 

1 1 



1 / 1 1_ 

p'l" - epb^ +iQE'- p'O - epb^ + zO " e' ~ e'^^ ^P? - ^ph +iO^ E'- pf - epb, +iOJ' 



(23) 



we obtain for the irreducible part 



r'^Y""''^ ^ — (l dE' i dE c/^^f''' (E' , E) 



J- J) dE' J) dE ^ ^ 
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1 



E ^ pf ~ en{l - iQ) 



(P6in|/(p'0-p;)|aia2)}. (24) 



The exp ression in the curly braces of equation (|24|l is a regular function of E or E' when E sa E^a"^ and i?' w i?^"-* 
[13: 



(see [lOj for details). Calculating the residues and taking into account the identity 

1 1 



27rVa:: + iO -x + iO 

we find 



5{x) , (25) 



I](-l)''{ {Pbi\ea^A'}*\n)— {nPb2\I{epb, - ea,)\aia2) 



(la.irrcd) 

p 11 ^ fcPbg ^ fcn 



^ (P62|ea^Ay*|n)^3^ (P5in|/(epb, - Jjaias) } . (26) 



A similar calculation of the irreducible part of the diagrams shown in Fig. 2b yields 



I](-l)''{ {PbiPb2\Iispb,-Sa,)\na2) _ ^ _ (n|ea^v4f |ai) 

V {PbiPb2\Iiepb,-ea,)\ain)^ {n\ea^A'i*\a2)} . (27) 



For the reducible part of the diagrams shown in Fig. 2a we have 

- hi/'^'i dEg^^\E\E) 



= / di?' / dE — — ^ \y{-i)p(-LY r dpidp\ 
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X {Pbi\ea^A>}*\n){nPb2\I{p'2° - P§)|aia2) 
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E - ^Pi - ^Pbi + *0 E- pf - e„ + iO 
1 1 



f L 



:{Ph2\ea^Af\n){PbMl{Pi-pl)\aia2)] 



(28) 



Calculating the residues at E' — E^^ and E = E^a^ and using the identity (|25|) . we obtain 



(la, red) 
^7j^ ,&;a 



■(£a,-pO+zO)2 



(P6i|ea^Af |n)(nP&2|/(epb, - P2)|aia2) 
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+ ? L27r7-oo"'''(£a,-p? + *0)2 
(P62 |ea^v4f |n) (P6in|/(ep6, - P?)|aia2)] } 



(29) 



We have assumed that the unperturbed states a and b are described by one-determinant wave functions Q and Q . It 
implies that, in equation H29|l . we have to consider {Pb2,n) = (01,02) or (02,01) in first term and {Pbi,n) = (01,02) 
or (02, oi) in second term. Therefore, the reducible part contributes only in the case when the states o and b have at 
least one common one-electron state. In what follows, we assume oi = bi and 02 ^ &2- We obtain 



/oc 
Ju{b2\ea,A^;\a2) 



(0iO2|/(cj)|0i02) (0201 |/(t^) 10102) 



{lU - i0)2 (t^ - Aa - ^0)2 

where = Ea^ — £a^ . A similar calculation of the reducible part of the diagrams shown in Fig. 2b gives 
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(31) 



where A^ = — £h^. The reducible contribution has to be considered together with second term in formula l|18|l . 
Taking into account that 
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(32) 



and a similar equation for the final state, one finds 

- i £^ dE' £ dE iE', E) £ dE ,« (^) + _L £ dE gil^ (E) 
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Summing (jSOI), ijSIJ, and we obtain for the total reducible contribution 



(33) 
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+ {b2bi\I{co)\bib2) 
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Here we have employed the symmetry property of the photon propagator: I{uj) = I[—uj). Using the identity 

1 1 2tt d 



i dtu 



and integrating by parts, we find 



(34) 



(35) 



(36) 



where /'(A) 



dljuj) 
duj 



uj=A 



and it is imphed that ai — hi. The total expression for t^^^j^.^ (in the case oi =bi) is given 



by the sum of equations ((23), l|2ZI), and 



^ , (la,irrcd) (16,iiTcd) (l.rcd) 



(37) 



In addition to the interelectronic-interaction correction derived above, we must take into account the contribution 
originating from changing the photon energy in the zeroth order transition probability due to the interelectronic- 
interaction correction to the energies of the bound states a and h. It follows that the total interelectronic-interaction 
correction to the transition probability of first order in 1/Z is given by 



(1) 



7/ 



(0) 



7/ 1^;*^ 



dW. 



(0) 



7/ ^^j*^ 



(38) 



where and E^'^ are the energies of the bound states a, h with and without the interelectronic-interaction 

correction, respectively. 



III. NUMERICAL RESULTS AND DISCUSSION 

To evaluate the one-electron transition matrix elements, the explicit formulas given in have been used. Infinite 
summations over the electron spectrum in equations (|26|l and 1)27(1 have been performed by using the finite basis set 
method. Basis functions have been constructed from B-splines by employing the procedure proposed in All the 
calculations have been carried out for the homogeneously charged sphere model of the nuclear charge distribution. 
The values for the nuclear radii were taken from (12] . 

In Tables HI Inl and lllll we present our numerical results for the decay rates of the magnetic transitions 2^ Si l^-S'o, 
2^P2 l^S'o, and 3'^Si — > 2'^ 5*1, respectively. The values presented in the upper and lower parts of the tables have 
been obtained in the Feynman and Coulomb gauges for the photon propagator, respectively. The transition energies 
used in the calculation were taken from [2, Q . The contribution due to the frequency dependence of the photon 
propagator (AM^ftoq) and the negative-continuum contribution {AW^+e-) are given in these tables as well. It can be 
seen from the tables that the total values of the transition probabilities in the different gauges coincide with each 
other. 

As one can see from Tables III and Ull for the decays with AS* ^ 0, the frequency-dependent correction is of the 
same and even larger magnitude than the negative-continuum contribution. However, this is not the case for the 
3'^S'i 2^Si transition, where the correction AWfrcq is small compared to the AW^+f,- term. The behavior of the 
negative-continuum correction as a function of the nuclear charge number Z agrees well with the scaling ratio of the 
negative- to positive-energy contributions found in Q for all the transitions under consideration. 

In Tables Ifvl and IVl we compare our results with the previous calculations 2, 3, 4] that partially include the 
and higher order terms but do not account for the frequency-dependent contribution. In Table Hvl the experimental 
data for the most precisely measured transition 2'' 5*1 I^Sq are also presented. In the last cohuim of this tabic our 
results are combined with the radiative corrections that are beyond the ones already included in the transition energy. 
These corrections were recently evaluated in J^J for the 2si/2 — *■ lsi/2 transition in hydrogenic ions for Z > 50. Since 
we consider high-Z two-electron ions, we can assume that the one-electron (hydrogenlike) approximation is sufficient 
to evaluate the related correction in He-like ions. We have extrapolated these data for Z < 50 and interpolated 



for Z = 54. The uncertainties due to the extrapolation of the radiative corrections and uncalculated 1/Z^ and 
higher order terms are indicated in parentheses. In Table El the comparison with the RMBPT calculations 4] is 
presented for the transitions 2^P2 ^^Sq and 3'^ Si 2^ Si. The uncertainties due to uncalculated radiative and 
higher order interelectronic-interaction corrections are also indicated. From Tables I and IV, it can be seen that the 
frequency-dependent contribution is smaller than the current experimental accuracy. 

In summary, we have presented a systematic quantum electrodynamic theory for the interelectronic-interaction 
corrections of first order in 1/Z to the transition probabilities in heliumlike ions. The numerical evaluation of these 
corrections to the magnetic transition probabilities has been performed and the equivalence of the Feynman and 
Coulomb gauges has been demonstrated. The results of the calculations performed have been compared with previous 
RMBPT calculations and with experiment. 
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TABLE I: The decay rates of the magnetic dipole transition — > l^So in units s^^. The negative-continuum contribution 
AWg+g- and the frequency-dependent correction AWfrcq are expressed in % with respect to the main term W. Wtot is the 
total decay rate value. The values presented in the upper part of the table were calculated in the Feynman gauge, whereas the 
results presented in the lower part were obtained using the Coulomb gauge. 

Z W AW^+^- AWfreq Wtot 

30 8.9994 X 10* -0.043% -0.029% 8.9929 x 10* 

50 1.7303 X 10" -0.08% -0.042% 1.7282 x lO" 

70 5.9872 X 10^^ -0.132% -0.045% 5.9766 x 10^^ 

90 9.4551 X 10^-^ -0.205% -0.036% 9.4323 x lO" 

30 9.0012 X 10* -0.05% -0.042% 8.9929 x 10* 

50 1.7308 X 10" -0.09% -0.062% 1.7282 x lO" 

70 5.9896 X 10^^ -0.145% -0.073% 5.9766 x 10^^ 

90 9.4596 X lO" -0.218% -0.070% 9.4323 x lO" 



TABLE II: The decay rates of the magnetic quadrupole transition 2^P2 — > 1^5o in units s~^. The negative-continuum 
contribution AW^+^- and the frequency-dependent correction AWfreq are expressed in % with respect to the main term W. 
Wtot is the total decay rate value. The values presented in the upper part of the table were calculated in the Feynman gauge, 
whereas the results presented in the lower part were obtained using the Coulomb gauge. 

z ir Air, I, AUfa.i Uu.t 

30 2.1047 x 10"' -0.0001% 0.021% 2.1052 x 10^° 
50 1.3654 X 10^^ -0.001% 0.038% 1.3660 x 10^^ 
70 2.1480 X 10" -0.005% 0.063% 2.1493 x lO" 
90 1.7231 X 10^^ -0.017% 0.097% 1.7245 x 10^'* 
30 2.1051 X 10^° -0.0001% 0.001% 2.1052 x 10^° 
50 1.3659 X 10^^ -0.001% 0.005% 1.3660 x 10^^ 
70 2.1491 X 10" -0.005% 0.014% 2.1493 x lO" 
90 1.7242 X 10" -0.017% 0.033% 1.7245 x lO" 
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TABLE III: The decay rates of the magnetic dipole transition 3^51 — > 2^ Si in units s~^. The negative-continuum contribution 
AWg+g- and the frequency-dependent correction AWfreq are expressed in % with respect to the main term W. Wtot is the 
total decay rate value. The values presented in the upper part of the table were calculated in the Feynman gauge, whereas the 
results presented in the lower part were obtained using the Coulomb gauge. 

Z W AW,+,- AW/frcq Wtot 

30 6.1245 X 10^ 3.867% 0.022% 6.3626 x 10=^ 

50 1.3019 X 10* 2.204% 0.034% 1.3311 x 10* 

70 4.9886 X lO'^ 1.488% 0.046% 5.0651 x lO'' 

90 9.0496 X 10^° 1.055% 0.059% 9.1503 x 10^° 

30 6.1273 X 10^ 3.837% 0.004% 6.3626 x 10^ 

50 1.3029 X 10* 2.158% 0.006% 1.3311 x 10* 

70 4.9936 X lO'' 1.428% 0.005% 5.0651 x lO'' 

90 9.0610 X 10"' 0.984% 0.002% 9.1503 x 10^° 



TABLE IV: The decay rate (s~^) of the transition 2^Si So calculated in this work is compared to the previous calculations 
and experiment. The experimental values and their error bars are given in second and fourth columns, respectively. In the last 
column the sum of our results and the QED corrections obtained in are presented. In parentheses the uncertainties of the 
present calculations are indicated. Relative differences are calculated using experimental results as a reference. 



z 


Exp. 


Ref. 


Prec. 


RMBPT \2] MCDF [3] 


Present Present -l-QED 


23 


5.917 X 10^ 


[14] 


4.1% 


-0.1% 


-0.4% 


0.1% 


0(6)% 


26 


2.083 X 10* 


[14] 


12.5% 


-0.4% 


-0.7% 


-0.3% 


-0 4(5)% 


35 


4.462 X 10^ 


[15] 


3.2% 


-2.3% 


-2.5% 


-2.1% 


-2.3(4)% 


36 


5.848 X 10^ 




1.3% 


-0.4% 


-0.6% 


-0.3% 


-0.5(4)% 


41 


2.200 X 10^° 




0.4% 


0.8% 


0.6% 


0.9% 


0.7(4)% 


47 


8.969 X 10^° 




1.8% 


1.3% 


1.2% 


1.4% 


1.1(2)% 


54 


3.915 X 10" 


m 


3.0% 


-1.8% 




-1.6% 


-2.1(2)% 



TABLE V: The decay rates (s ^) of the transitions 2^P2 l^So and 3^51 — » 2^51 obtained in this work are compared to the 
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FIG. 3: One-photon exchange diagram. 



